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Abstract 

The electromagnetic propagator associated with the localized collective mag- 
netoplasmon excitations in a semiconductor superlattice with broken transla- 
tional symmetry, is calculated analytically within linear response theory. We 
discuss the properties of these collective excitations in both radiative and non- 
radiative regimes of the electromagnetic spectra. We find that low frequency 
retarded modes arise when the surface density of carriers % at the symmetry 
breaking layer is lower than the density r] at the remaining layers. Otherwise 
(Vd > ? 7o) a doublet of localized, high-frequency magnetoplasmon-like modes 
occurs. The corresponding dispersion law and power spectrum of these modes 
are shown. 
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I. INTRODUCTION 



Much work has been devoted to the study of the properties of low- dimensional het- 
erostructures [Q. Many of the advances are due to the emergence of new crystal growth 
techniques, such as molecular-beam epitaxy (MBE) and metal organic chemical vapour de- 
position (MOCVD). These experimental techniques have made possible the fabrication of 
layered materials with sharp, high-quality interfaces and with dimensions comparable to the 
electron mean free path and the de Broglie wavelength ||[§. There are many experimental 
and theoretical reports concerning the collective excitations that can propagate in such struc- 
tures, including magnons, plasmons, magnetoplasmons and polaritons [4-17]. Experimental 
set-ups for detecting some of these collective oscillatory modes have been reported recently 
|T6| , p!7[| . Most of the theoretical work in this field is concerned with infinite and semi-infinite 
superlattices (for a review see [fR|). For a periodic array of quantum wells, a plasmon band 
arises because of the translational symmetry and the long-range electromagnetic coupling 
of the two-dimensional (2D) electron sheets [^P,|i~9f. The spectrum of the collective modes 
becomes richer when an external static magnetic field is applied along the superlattice axis. 
In particular, low-frequency undamped helicons ||T0| , |20| -[22| and waveguide-like modes |23| 



can be detected. The case of a semi-infinite semiconductor superlattice was considered by 
Giuliani and Quinn ||IHJl2|j2"4l . They found a polariton-type plasmon mode (Giuliani-Quinn 
polariton) which exists only for values of the in-plane wave-vector greater than a certain 
critical value. This critical value is determined by the differences between the dielectric 
permittivities of the bulk media and the semi-infinite superlattice. The extension of these 
works to the case of applied magnetic fields was discussed by Kushwaha p3| . Recently much 
attention has been paid to the description of the local modes related with the removal of 
translational symmetry, as for a superlattice with a defect layer [|l3| , p6f or a superlattice with 
a quasiperiodic region following Fibonacci |27) or Cantor |28| sequences. The effects of an 
applied external magnetic field have also been considered using Green's function techniques 

Studies of magnetoplasmons in layered heterostructures have mainly focussed on the dis- 
persion law of collective excitations, however the power spectrum of such excitations has 
not been discussed to our knowledge. In the present work we calculate the retarded pho- 
ton Green's tensor for the magnetoplasmon modes arising in superlattices with broken 
translational symmetry. We also discuss the dispersion law (by taking into account retarda- 
tion effects and the non-homogeneous propagation of these modes) and the power spectra 
of low-frequency collective excitations. 

The outline of the paper is as follows: In Section II we describe the system considered 
here and obtain analytical expressions for the Green's functions associated with the localized 
collective oscillations. The corresponding discussion, which is centered on the properties of 
the collective excitations in both non-radiative and radiative regimes of the electromagnetic 
spectra, is given in Section III. Analytical and numerical treatments are presented for the 
cases of non-retarded (Sec. Ill A) and homogeneous (Sec. Ill B) modes of localized mag- 
netoplasmon oscillations. The behavior of the retarded non-homogeneous localized modes is 
discussed in Sec. Ill C. Conclusions are given in Sec. IV. 
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II. FORMALISM 



We consider an infinite one-dimensional array of quantum wells forming a Type-I super- 
lattice with periodicity d. The z— direction is taken to lie along the superlattice axis. We 
assume that all the quantum wells are uniformly doped with a surface density of carriers 
r] , except the one centered at z = which is doped with a surface density r) d . An external 
static magnetic field B is applied along the z— direction. We also assume that the distance d 
between wells is so large that we can neglect wave-function overlap between wells. Therefore 
the quantum wells are coupled only by the electromagnetic interaction associated with the 
dynamics of the electron system. 

The dynamic magnetoconductivity tensor corresponding to the n'th layer is given by 



[n) 



<4 + 



nO) 



where (afj) are the components of the 2D dynamic magnetoconductivity tensor associated 
with the electronic planes with surface density rj The photon Green's tensor Dj k 

corresponding to the dynamics of the electron system satisfies the following set of differential 
equations: 



d 2 



dxidxj 



+ V 2 8i^j D jk (ou;r n ,r' v z,z') =4^^5(7--^)- 
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where co is the frequency of the allowed collective excitations, and r' n are vectors in the 
xy— plane, r = (r lh z) and A(r] d ) — — 1. The superscript (n) on the tensor D jk denotes the 
layer under consideration. Given the homogeneity of the system along the electron layers, 
we introduce the 2D Fourier transformation 



D jk (u;r,r') = -^-^ J d 2 k }l d jk (uj, fe„, z, z')e 



ik r (r,, - r' ) 
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where fey is the in-plane wave vector with components (k x , k y ). The isotropy in the xy— plane 
allows us to perform a coordinate rotation which aligns the x axis with the direction of k n . 
This is achieved by the action of the matrix 
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Equation (2) hence takes the form 
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where — Sj m d mn S n ^ . In this matrix relation, we introduce the notation gj k = 
g jk {uj, k n ;z, z'), g$ = g jk{n) = g jk (w, kf,z = nd, z') and k 2 = k 2 - *£e. 

Without loss of generality, we consider the solutions of Eq. (5) in the region jd < z < 
(j + l)d, with j 7^ and < z' < d. The solutions can be represented as follows: 
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satisfying the following boundary conditions at the n'th layer 
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with 
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To solve for the local modes in the superlattice considered here, we assume a solution 

sS(n) = »S(o) exp(-a|n|d), (11) 

which decays exponentially as n — > ±oo, i.e. far from the electron layer with density r/d 
(z = 0). Equation (11) shows that the attenuation of the fields along the superlattices axis 
is characterized by the parameter a. The condition for the existence of non-trivial solutions 
of Eq. (5) in the regions considered above, yields 

\ 1xx S{k, a) + l) (^- lyy S{K, a) + l) = -1^S 2 (k, a), (12) 



where 
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is the superlattice structure factor. Equation (12) is expressed in a general form in terms of 
the magnetoconductivity tensor er^-. We employ the local approximation for Oif 
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is the cyclotron frequency of the two-dimensional carriers. In this approxi- 
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mation, Eq. (12) takes the form 
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-f- is the three-dimensional plasmon frequency. Equation (15) would be 



where u v — , 
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equivalent to that obtained in Refs. |9||n| if ^ a were to be replaced by a quasi wavenum- 
ber. When each quantum well is doped with the same surface density, a magnetoplasma 
band arises which is related to the periodicity of the system along the superlattice direction. 
However if rjd ^ i] (i.e. if a is a complex number), then | cosh(ad)\ > 1: localized magne- 
toplasmon modes associated with the broken translational symmetry appear on the outside 
of the bulk magnetoplasmon band. 

Following a similar procedure to that described in Ref. | 59fl , it is straightforward to 
show that the diagonal components of the Green's tensor in the region — d < z < d, with 
< z' < d, have the following form: 
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In these expressions we have introduced the shorthand notation: 
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Notice that the poles of Eqs. (16-18) together with Eq. (15) give 
the allowed frequencies of the modes which are localized at the defect layer. They can be 
found analytically in the non-retarded limit and for homogeneous oscillations in the retarded 
region. These cases, together with retarded non-homogeneous oscillations, are discussed in 
the next section. 



III. RESULTS AND DISCUSSION 



A. Non-Retarded Localized Modes 



ck„ 



In the non-retarded region of the spectrum uj « and the components of the Green's 
tensor have poles at: 
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Here A (77) = 6(77 — 1) — 6(1 — 77), where Q(x) is the Heaviside step function. In order to have 
a solution describing the local magnetoplasma mode, we need X = cosh(ad) > 1 or X < 1. 
In the non-retarded limit, expression (35) does not depend on the external magnetic field and 
a solution of X exists for all values of the in-plane wave- vector k n . This behavior contrasts 
with that of the Giuliani-Quinn surface polariton |T]J which can not be excited in the long 
wavelength region of the spectrum. The above general features are illustrated in Figs. 1 and 
2, where we have plotted the dispersion curves u = co(k n ) of the localized magnetoplasma 



mode for different values of rj. The regions between the dashed lines correspond to the 
allowed frequencies of the magnetoplasmons in the ideal system. It can be seen that the 
local frequency increases its value with k n . When k n d ^> 1, the frequency of the local mode 
approaches the frequency of the magnetoplasmon oscillations of a single quantum well with 
surface carrier density This result is reasonable since the coupling between quantum 
wells is weak in the k^d >> 1 limit. Figure 1 corresponds to values of r\ < 1 : in this case 
X < — 1 and the local frequencies are lower than those corresponding to the edge of the 
Brillouin minizone (X = — 1). The lowest frequency in this case corresponds to ui c . In the 
strong coupling limit (k n d <C 1), the local mode has a frequency whose dispersion law can 
be expressed as 

u = uj c + ^fr 1 (l-X){k,d)\ (36) 



ed ' V 2 
where 

^4(H (M)2 -(H' (37) 

with a H = In other words, the behavior of the local low- frequency, non-retarded modes 
in the strong coupling limit is governed by the Hall conductance of the 2D electron gas a H . 
We might conclude from this analysis that the frequency of this mode is always higher than 
the cyclotron frequency of the two-dimensional carriers. We will see however that accounting 
for retardation effects can alter this result drastically. 

Figure 2 corresponds to rj > 1. Then X > 1 and the local frequencies are higher than 
those corresponding to the center of the Brillouin minizone (X = 1). In the strong coupling 
limit we have 

(38) 



2 2,2 

u =u c +u p 



i + Q + f^- 2 )) (M) 



where 

e ad = l + ^-l)(M) 2 - (39) 

We see t hat the frequency of this local mode approaches the three-dimensional value 
Up + u 2 c as fc|| approaches zero. 

B. Homogeneous Magnetoplasmon Oscillations 



In the above discussion we have neglected the role of retardation effects. Now we consider 
these effects for the k^ = modes, which corresponds to homogeneous (along the planes) 
oscillations of the electron system. In this case 



cot (^] + Afo) {cot 2 )- v ( V - 2)\ , (10) 




u = ±u Jc + — 

e ad ={l-^os{^)+A( V ){[l--J co^[^) + [~ z -l)} . (41) 

In contrast with previous non-retarded results, the magnitude e ad in the present case shows 
a dependence on the external magnetic field B . This follows from the fact that the modes 
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described by Eq. (41) owe their existence to retardation effects. In this sense the origin of 
these local modes differs markedly from the origin of the non-retarded ones. In the particular 
case when |1 — rj\ <C 1 we find localized modes with frequencies close to the center or the 
edge of the Brillouin mini- zone (X = ±1 respectively). We obtain for these frequencies that 

u = u - A( V ) (1 - ,f ^sin , (42) 

where ui Q is the frequency at which the local mode enters the bulk band of the waveguide-like 



magnetoplasmon oscillations which exist in the ideal system |23| . For u> c > the frequency 
uj a satisfies the inequality ui < uj c . This means, that the local modes have frequencies lower 
than the electron cyclotron frequency. The above features are illustrated in Figs. 3 and 4, 
where we have plotted the dependences ui = u)(rj) for different values of the applied external 
magnetic field. 

Figure 3 shows the dependence ui = u>{rj) for the case i] d < rj Q . When r] d tends to zero, the 
local frequency approaches the cyclotron frequency. This means that the absence of carriers 
at the defect layer leads to the appearance of a local mode with a frequency equal to the 
cyclotron frequency of the carriers at the remaining layers. We also see that for a fixed value 
of the surface carrier density rjd at the defect layer, the local frequency increases with the 
applied magnetic field. For a given value of the external applied magnetic field Bo, however, 
to decreases with increasing rjd. In the case rjd > f] the frequencies corresponding to localized 
modes appear only in the high frequency retarded region (i.e. uj ^> uj c ). The analysis shows 
that there exist a large number of doublets of local modes. Figure 4 displays the dependence 
uj = uj{rj) for the first doublet. For a given value of rjd we found that the frequency of the 
lowest (highest) component of this doublet decreases (increases) with applied magnetic field. 
For this reason the separation between the components of the doublet increases considerably 
with magnetic field. 



C. General Dispersion Relations and Power Spectra 

We now focuss on the description of the localized modes in the retarded region 
(i.e. uj > Figure 5(a) shows the behavior of the frequency uj (in units of ui p ) as a 

function of in-plane wave number (in units of k) f° r V = 0.5. The dashed region corre- 
sponds to the bulk helicon band. We see that in the forbidden region, a local helicon mode 
appears which is localized in the vicinity of the bulk band edge. With increasing in-plane 
wave number k n , the frequency of the local modes rises until it reaches the electron cyclotron 
frequency; here it loses its local-mode character since its frequency is now indistinguishable 
from the frequencies describing the layer dynamics in the ideal system. Also we note that if 
r] — > 1, the dispersion curves aproach the edge of the bulk helicon band; if on the other hand 
r] —>■ 0, it can be seen from the poles of the Green's functions (Sec. II) that the frequency 
corresponding to the dynamics of the defect layer tends to ui c for all values of k n . In Fig. 5(b) 
the first doublet of localized modes corresponding to rj = 1.5 is plotted. The dashed area 
corresponds to the bulk magnetoplasma band. We see that the lowest mode of the doublet 
slowly increases its frequency with k n until it reaches the bulk continuum. The other mode 
of the doublet increases its frequency more rapidly than the lowest mode, and tends to the 
dispersion relation u = with increasing fcy. 



The expressions for the diagonal Fourier components of the photon propagator enable us 
to obtain the peaks corresponding to power spectra of such modes; i.e. the imaginary part of 
the sum over gu given by Im^ 9u where i = x,y, z. In Fig. 6(a) we show the power spectra 
of the lowest frequency modes for different values of fey, in the regime where the density 
of carriers at the defect layer satisfies % < r] . Only one peak exists which corresponds to 
frequencies lower than the cyclotron frequency uj c . With increasing in-plane wave vector , 
the frequency at the maxima of such peaks shifts to the right until it reaches u c . These 
features are also illustrated in Fig. 5(a). In addition, the frequency corresponding to the 
peaks in Fig. 6(a) increases when the density of carriers at the defect layer diminishes. 

For the case rjd > rj , there are no peaks with frequencies satisfying u < uj c . In this case 
the peaks corresponding to localized modes appear in the high frequency retarded region of 
spectrum. The analysis shows that there exist a large number of doublets of local modes. 
We limit ourselves to discussion of the features corresponding to the first doublet. In Fig. 
6(b) we show the power spectra for different values of k n . In the homogeneous case the 
separation between the components of the doublet is small in comparison with uj p and it 
increases with increasing in-plane wave vector. These results agree with those in Fig. 5(b). 



IV. CONCLUSIONS 



We have used the formalism of linear response theory in order to consider the localized 
modes in a Type-I superlattice with broken translational symmetry, in the presence of an 
external static magnetic field. We found that a low frequency retarded mode arises only when 
the surface carrier density rjd is lower than the corresponding density rj Q at the remaining 
layers. This mode has frequencies greater than those corresponding to the edge of the 
helicon band, but is lower than the cyclotron frequency. In contrast for rj d > rj , we found 
a doublet of localized high-frequency magnetoplasmon-type modes. These results can be 
extended to consider a distribution of defect layers in a Type-I superlattice. In this case it 
is necessary to take into account the techniques developed for the description of localized 



states in disordered systems To our knowledge there are so far no experimental results 
on superlattices with a single defect layer, however we believe that the localized collective 
excitations described here can indeed be detected with currently available technology using 
Raman spectroscopy experiments: similar experiments have been reported recently for the 
case of GaAs/Al x Gai_ x As double quantum well structures jlSfl , where acoustic and optic 
plasmon modes of this electron bilayer system were successfully observed using electronic 
Raman scattering. We hope that the results reported in this paper will stimulate further 
experimental work on this topic. 

J.H.R. thanks the financial support of COLCIENCIAS. 
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Figure Captions 



FIGURE 1. Dispersion curves of the local non- retarded magnetoplasmon modes for 
different values of the surface carrier density at the defect layer: 77 = ^ = 0.1, 0.5, and 
0.9. The regions between the dashed lines correspond to the allowed frequencies of the 
magnetoplasmons in the ideal system. u c = 0.014^ and = 2.3. 

FIGURE 2. Same as in Fig. 1 but for 77 = 1.5, 77 = 3.0, and 5.0. 

FIGURE 3. Plot of the homogeneous (k n = 0) low-lying local magnetoplasmon fre- 
quencies as a function of the surface carrier density at the defect layer for 77^ < rj and for 
different values of the applied magnetic field: u c = 0.014u;p (Bq = IT), 0.14a;p (Bo = 10T), 
and 0Alu p (B = 30T). = 2.3. 

FIGURE 4. Same as in Fig. 3 but for r] d > r] . 

FIGURE 5. (a) Dispersion curves of the localized magnetoplasma modes (solid lines) for 
retarded non-homogeneous oscillations (u > ^0 . 77 = 0.5, u c = 0.014^ and Upd ^ = 2.3. 
The dashed region corresponds to the bulk helicon band, (b) Same as in Fig. (a) but for 
77 = 1.5. The dashed area corresponds to the bulk magnetoplasmon band. 

FIGURE 6. (a) Power spectra (in arbitrary units) of low-frequency localized magneto- 
plasmon modes as a function of the excitation frequency and the in-plane wave number in 
the cases k^d = 0.0, 0.015 and 0.025. Input parameters as in Fig. 5(a). (b) Same as in Fig. 
(a) but for k^d = 0.0, 1.0 and 1.5. Input parameters as in Fig. 5(b). 



1 1 




Figure 1 

J.H. Reina et al. 




Figure 2 

J.H. Reina et al. 




Figure 3 

J.H. Reina et al. 




Figure 4 

J.H. Reina et al. 




Figure 5(a) 

J.H. Reina et al. 




Figure 5(b) 

J.H. Reina et al. 




Figure 6(a) 

J.H. Reina et al. 




Figure 6(b) 

J.H. Reina et al. 



